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We recently questioned the treatment of a dust particle as a perfect absorber for electrons and ions and proposed 
a surface model for the charge of a dust particle in a quiescent plasma which combines the microscopic physics 
at the grain boundary (sticking into and desorption from external surface states) with the macrophysics of the 
discharge (plasma collection fluxes). Within this model the charge and partial screening of the particle can be 
calculated without relying on the condition that the total electron collection flux balances on the grain surface 
the total ion collection flux. Grain charges obtained from our approach compared favorably with experimental 
data. The purpose of this paper is to describe our model in more detail, in particular, the hypotheses on which 
it is built, contrast it with the standard charging models based on flux balancing on the grain surface, and to 
analyze additional experimental data. 
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1 Introduction 

Dust particles immersed in an ionized gas acquire charges whose signs and magnitudes strongly affect the prop- 
erties of the host plasma as well as the properties of the ensemble of dust particles itself. The various crystalline 
and liquid phases of dust particles found in specifically designed laboratory experiments are perhaps the most 
impressive manifestations of this effect |[Tl|2l. That dust particles and host plasma are strongly coupled has 
been however known for a long time. Not only from astrophysical plasma environments 13] ID but also from 
processing discharges, where the grains are of course not the intentional constituents studied for their own sake 
but contaminants which need to be controlled because they are detrimental to the overall performance of the 
discharge IS]. 

The charge of a dust particle is an important parameter for the quantitative description of dusty plasmas. It has 
thus been measured in a number of experiments ||6l|7][8]|9][l0l[TTl[l2l. Considering the dust particle as a floating 
electric probe, the interpretation of the experiments is based on two main assumptions: (i) The grain is a perfect 
absorber for electrons and ions, that is, every electron and ion hitting the surface of the grain is absorbed and (ii) 
the quasi- stationary charge of the dust particle is the one which balances on the grain surface the total electron 
collection flux with the total ion collection flux |fT3l [T4l [T5l [161 [TTI [TSl [191 . The charge enters here through 
the floating potential of the dust particle which, for a spherical particle, can be immediately translated into a 
charge. In most cases, however, the charges obtained from this approach are not in agreement with experimental 
data. Usually, the approximations for the fluxes are blamed for the disagreement |[T6l [TSl [191 [20l [211 [22l [23l . 
We suspect, however, that the flux balance condition, as it is currently used, is incomplete because it ignores 
important processes on the grain surface which also affect the grain charge. 

First indications that this could be indeed the case came from a phenomenological analysis of experimental 
data based on rate equations for the electron and ion surface densities of a dust particle 1241 . To provide a 
more realistic model for the grain surface, the rate equations contained electron and ion sticking probabilities, 
electron and ion desorption times, and a electron-ion recombination probability. With an appropriate choice of 
these parameters as well as the grain temperature the charges obtained from the quasi-stationarity of the surface 
densities were in better agreement with experimental data then the ones deduced from the orbital motion limited 
flux balance condition. 
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Fig. 1 Left panel: Generic illustration of elementary surface processes leading to the formation of surface charges at an 
inert plasma boundary where surface modifications other than charge accumulation do not occur. Right panel: Surface model 
proposed in 1251 for the region around a dust particle with radius R. At quasi-stationarity, the surface charges (Te,, bound at 
re — R and r; > re, respectively, balance the collection fluxes Se,ije^^^"^^ with the respective desorption fluxes r^lae,i, 
where s^^i and r^^i denote, respectively, sticking coefficients and desorption times. 



In a semi-microscopic approach we subsequently pointed out that the charging of a dust particle can be inter- 
preted as a physisorption process in the potential which builds up in the disturbed region around the grain ||25l . 
The potential is the sum of a short-range polarization potential which enforces the electric boundary conditions on 
the grain surface and, as soon as the grain collected some charge, a long-range Coulomb potential. Realizing that 
the analogy to physisorption suggests, on a microscopic scale, a spatial separation of bound electrons and ions, 
we calculated the charge of the grain and its partial screening by balancing, on two different effective surfaces, 
the electron collection flux with the electron desorption flux and the ion collection flux with the ion desorption 
flux. The charge of the dust particle thereby obtained is then given by the number of electrons quasi-bound in the 
polarization-induced short-range part of the particle potential whereas its partial screening is given by the ions 
quasi-trapped in the long-range Coulomb part of the particle potential. 

In this paper we discuss the surface model proposed in ||25l, in particular, the hypotheses on which it is built, 
in more detail and contrast the model with the standard approaches of calculating the particle charge from flux 
balancing on the grain surface. Thereby we identify certain issues which need to be resolved before a more 
refined microscopic theory of particle charging can be developed. Finally, to demonstrate that the surface model 
in its present form works we analyze further experimental data and show that the charges obtained from it are 
much better than the ones deduced from the standard approaches. 



2 Elementary surface processes 

The microphysics at the grain surface affecting the charge of the grain is schematically shown in the left panel of 
Fig[T] Electrons and ions are collected from the plasma with collection fluxes = Se.ijg'f'''"'', where Se,i are 
the sticking coefficients and j^^^"^'^ are the fluxes of plasma electrons and ions hitting the grain surface. Electrons 
and ions may thermally desorb from the grain surface with rates r~/, where Te,i are the desorption times. They 
may also move along the surface with mobilities /ig.i, which in turn may affect the probability an with which 
ions recombine with electrons on the surface. All these processes occur in a thin layer whose thickness is at 
most a few microns, that is, on a scale where the standard kinetic description of the gas discharge based on the 
Boltzmann-Poisson system to which the flux balance condition belongs breaks down. 

To put the surface model proposed in Ii25il into context and to identify the assumptions with respect to the 
surface properties which are usually made in the standard calculations of surface charges we cast the elementary 
surface processes into rate equations. 
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Specifically, we consider a spherical dust particle with radius R. The quasi-stationary charge of the grain is 
given by (we measure charge in units of — e) 

Zp = AnR^ap = AttR^ [ae - (Ti] , (1) 

with electron and ion surface densities, cte i, satisfying the quasi-stationary {d(Je,i/dt = 0) rate equations ||24| . 

= SejPi— - T- V, - anaeu, , (2) 

-plasma — 1 /o^ 

= SiJi - CTi~ anaeCTi . (3) 

where j^'f*'™'^, Sg i, r^^i, and have been introduced above. For simplicity we neglect diffusion along the 
surface. 

In order to derive the standard criterion invoked to determine the quasi-stationary grain charge, we now as- 
sume, in contrast to what we do in our model |25| (see below), that both electrons and ions reach the surface of 
the grain. In that case, both Eq. ^ and Eq. ^ are flux balances on the grain surface. At quasi-stationarity, the 
grain is charged to the floating potential. In energy units U = ZpC^ /R. Because the grain temperature kTp <C U 
the ion desorption rate ^ ~ 0. Equation Q reduces therefore to fie (Tj = s jf''*^™^ which transforms Eq. (|2]) 
into Seje^^^'^^^ = Sijf^'^^™'^ + T~^(jp provided cTp ~ ae which is usually the case. In the standard approach the 
grain surface is assumed to be a perfect absorber for electrons and ions. Thus, Si = s,. = 1 and = 0. The 
quasi-stationary charge Zp of the grain is then obtained from the condition 

jP'— (Zp) = jt'"'^{Zp) , (4) 

where we explicitly indicated the dependence of the plasma fluxes on the grain charge. 

Calculations of the grain charge differ primarily in the approximations made for the plasma fluxes j^'f For 
the repelled species, usually coUisionless electrons, the flux can be obtained from Poisson's equation and the col- 
lisionless Boltzmann equation, using trajectory tracing techniques based on Liouville's theorem and momentum 
and energy conservation lfT3l [T4l [iTl . The flux for the attracted species, usually collisional ions, is much harder 
to obtain. Unlike the electron flux, the ion flux depends not only on the field of the macroscopic body but also 
on scattering processes due to the surrounding plasma, which throughout we assume to be quiescent. For weak 
ion collisionalities the charge-exchange enhanced ion flux model proposed by Lampe and coworkers 120112111221 
is often used. Its validity has been however questioned by Tskhakaya and coworkers 1261 |27l . We come back to 
Lampe and coworkers approach below when we discuss representative results of our model. 

Irrespective of the approximations made for the plasma fluxes, the standard approaches of calculating surface 
charges are based on three assumptions about the surface physics: 

(i) Ions and electrons reach the surface, even on the microscopic scale, 

(ii) Se ~ Si ^ 1 or at least Se = s;, and 

(iii) = or at least t^^cTc ^ Sijf^^^"^^ ~ aRUef^i- 

Hence, electrons and ions hitting the surface of the grain are assumed to be completely absorbed. For low- 
temperature gas discharges, where average electron energies kT^ ~ 10 eV, this cannot be the case, however 
Permanent implantation of electrons with this energy is very unlikely. Instead, electrons impinging with a few 
eV on a surface are either reflected, inelastically scattered or, when the surface potential supports bound states, 
temporarily trapped in external surface states, with residence times depending on the inelastic coupling of the 
electrons to the elementary excitations of the surface and the bulk material. As far as ions approaching the dust 
particle are concerned they can gain at most the floating energy of a few eV . They are thus also not able to enter 
the grain and to get permanently stuck. Q 

Based on these considerations, the three assumptions concerning the grain surface made in the standard cal- 
culations of the grain charge seem to be rather unrealistic. We challenge therefore all three of them. 

First, electrons and ions should be bound in external surface states. Because of differences in the potential 
energy, mass, and size we expect the spatial extension of the electron and ion bound states, and thus the distance 
of electrons and ions from the grain surface, to be different. On the microscopic scale, electrons and ions should 
be spatially separated. 



' We do not account for ion neutralization due to electron capture within the disturbed zone of the grain. 
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Second, Sg = Si is quite unlikely. Usually, the sticking coefficient of heavy particles is determined by the 
coupling to vibrational excitations of the surface and the bulk material 1281 |29l . This coupling is very strong. 
Thus, if ions reach the surface, as it is conventionally assumed, they would efficiently dissipate energy. The 
sticking probability would be thus large, that is, close to unity. Light particles, like electrons, on the other hand, 
couple only weakly to surface and bulk vibrations. We would thus expect Sg ^ s^. In principle the coupling 
to other elementary excitations of the grain (plasmons, electron-hole pairs, ...) can compensate for the lack of 
coupling to lattice vibrations but how efficient this coupling really is is not at all obvious. In a recent exploratory 
calculation we found, for instance, that the electron sticking coefficient at metallic boundaries arising from the 
coupling to internal electron-hole pairs is also extremely small, at least when the Coulomb interaction between 
an internal and an external electron is realistically screened 1301 . 

Third, if ions and electrons are indeed spatially separated, the two rate equations (without the recombination 
term) should be in fact interpreted as flux balances on two different effective surfaces as shown in the right panel 
of Fig. [T] In that case, ajjCTifTe <C o-e.i/ ''"e,i and the surface charge Zp would be determined by balancing on the 
grain surface the electron desorption flux, r^^fje, with the electron collection flux, Se.?^'**^™^. The correspond- 
ing balance of ion fluxes, to be taken on an effective surface surrounding the grain, would then yield a partial 
screening charge Zi. Within this scenario, we would thus obtain 

Zp = 47rr2.(sr),-jP'--^(Zp), (5) 
= 47r,^(,sr),.jf^^'°% (6) 

with re — R and > r^, which are the equations proposed in ||25l . The sticking into and desorption from 
external surface states is here encoded in the products (sr)e i. For electrons the product depends only on material 
parameters but for ions it also depends on plasma parameters. 



3 Surface states 

The starting point of our surface model was a quantum mechanical investigation of the bound states of a negatively 
charged particle in a gas discharge. More specifically, we considered the static interaction between an electron 
(ion) with charge — e (+e) and a spherical particle with radius R, dielectric constant e, and charge Qp = —eZp. 

As mentioned before, the interaction potential contains a polarization-induced part, arising from the electric 
boundary conditions at the grain surface, and a Coulomb tail due to the particle's charge ||3T1 [32l . For both 
terms we adopted the simplest approximations. The Coulomb part is then the potential of a point charge whose 
magnitude is the charge of the grain and the polarization part is the classical image potential. More sophisticated 
treatments taking, for instance, the finite velocity of the approaching electron (ion) and the nonlocality of the 
polarization potential at short distances into account are possible but not needed at this stage of the investigation. 

Measuring distances from the grain surface in units of R and energies in units of U, the interaction energy at 
X = r/R — 1 > Xf), where Xb is a lower cut-off, below which the grain boundary cannot be described as a perfect 
surface anymore, reads 

Ve,^{x) = ±-_ - 

1 — ^/2.T electron 

-1/(1 + a;) ion ^' 

with^= (e-l)/2(e+l)Zp. 

The second line in Eq. (|7]i is an approximation which describes the relevant parts of the potential very well and 
yet permits an analytical calculation of the surface states. In Fig.|2]we plot Ve.i{x) for a melamine-formaldehyde 
(MF) particle (e ~ 8, i? = 4.7 ^m, and Zp = 6800) embedded in a helium discharge with plasma density = 
rij = 0.62 X lO'' cm"^, ion temperature fesT, = 0.04 eV, and electron temperature ksTe = 2.2 eV ll24ll33l . 
From the electron energy distribution, fe{E), we see that the discharge contains enough electrons which can 
overcome the Coulomb barrier of the particle which is the floating energy U. These electrons may get bound in 
the polarization-induced short-range part of the potential, well described by the approximate expression, provided 
they can get rid of their kinetic energy. Ions, on the other hand, being cold (see fi{E) in Fig. |2]i and having a 
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Fig. 2 Left panel: Potential energy for an electron (ion) in the field of a MF particle (7? = 4.7 /xm, Z = 6800) I24II33I and 
representative probability distributions, |ii(a::)p, shifted to the binding energy and maxima normalized to one. Dashed lines 
denote the potentials used in the Schrodinger equations. Note, the finite ion radius pp^^ ~ A forces the ion wavefunctions to 
vanish at a; < il''*^ ~ 10"*. Right panel: Bulk energy distribution functions for the discharge hosting the particle I24II33I . 



finite radius p^'^'^/i? = xl*^*^ ~ 10~*, cannot explore the polarization potential at short distances 0. For them, 
the long-range Coulomb tail is most relevant, which is again well described by the approximate expression. 

Writing for the electron eigenvalue e*^ = 1 — Q;e^/4/c^ with ae — {e — l)i?/4(e + \)aB and for the ion 
eigenvalue = ~ai/2k'^ with = niiUZp/meaB, where as is the Bohr radius and me and are the 
electron and ion mass, respectively, the radial Schrodinger equations for the approximate potentials read 



e.i 
"fc2" 



l{l + l) 



(1 







(8) 



with Ve{x) = 2ae/x and Vi{x) = 2ai/{l + x). For bound states, the wavefunctions have to vanish for x ^ oo. 
The boundary condition at Xf, depends on the potential for x < x^, that is, on the surface barrier (which is 
different for electrons and ions). Matching the solutions for x < Xb and a; > Xf, at .t = x^ leads to a secular 
equation for k. For our purpose, it is sufficient to take the simplest model for the barrier: Ve,i {x < Xf,) = oo with 
Xb — for electrons and Xb = xf^^'^ for ions. 

The electron Schrodinger equation is then equivalent to the Schrodinger equation for the hydrogen atom and 
k is an integer n. Because (for bound electrons) a; ^ 1 and ^ 1, the centrifugal term is negligible. Hence, 
we consider only states with I = 0. The eigenvalues are then = 1 — ae^/4n^ and the wavefunctions read 



,,,ix) ~ ^;„,o(^-) = z-exp(-z72)(-)"-i(n - ml'^z) 



(9) 



with z = 2aex/n and il,^-' (z) associated Laguerre polynomials. 

The probability densities \u^ o(^)l^ '^^e first three electron surface states are plotted in Fig.|2] As can be 
seen, they are very close to the surface. In physical units, the electron surface states are at most a few Angstroms 
away from the grain boundary. At these distances, the spatial variation of Ve{x) is comparable to the de-Broglie 
wavelength of electrons approaching the particle. More specifically, for ksT^ = 2.2 eV, Ag^/i? ~ IK/K'I — 
10~"'. Hence, the trapping of electrons at the surface of the particle is a quantum-mechanical effect beyond the 
Boltzmann-Poisson description of the plasma-grain interaction. 



We treat the ion as a structureless rigid sphere. 
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The solutions of the ion Schrodinger equation are Whittaker functions, ul i{x) = Wk.i+i/2{x) with x = 
2(3;, (1 + x)/k and k determined from ((^cf ^'^) ~ 0. However, since ^ 1 and fc ^ 1, it is very hard to work 
directly with Wk,i+i/2{x). It is easier to use the method of comparison equations 13411 and to construct uniform 
approximations for u], ; (x) with the radial Schrodinger equation for the hydrogen atom as a comparison equation. 
The method can be applied for any /. Here we give only the result for I ~ 0: 

■"fc,o(-^) ~ Vn^oiz) /\/dz/dx (10) 

with Vnfl{z) defined in Eq. (|9]l and z = 2aiz{x) / n. The mappings z{x) and k{n) can be constructed from the 
phase-integrals of the two ion Schrodinger equations ll34l . 

In Fig.|2]we show |ufe.o(^)P for '^(1)' ^(5000) and fc(lOOOO). Note, even the fc(lOOOO) state is basically at the 
bottom of the potential. This is a consequence of ^ 1 which leads to a continuum of bound states below the 
ion ionization threshold at e = 0. We also note that o(^)l^ peaks for n':$' 1 just below the turning point. 

Hence, except for the lowest ion surface states, which we expect to be of little importance, ions are essentially 
trapped in classical orbits deep in the disturbed region of the grain. This will be also the case for / > 0. That ions 
behave classically is not unexpected because for fesT; = 0.04 eV their de-Broglie wavelength is much smaller 
then the scale on which the potential varies for x > lO^'^: Xf^/R ~ 10^^ ^ IK/V^'I — 1- Thus, the interaction 
between ions and the dust particle is classical and can be analyzed with Boltzmann-Poisson equations. 

Nevertheless it is also possible to describe this interaction quantum-mechanically. We anticipate even a 
quantum-mechanical approach, based on the method of comparison equations, which is an asymptotic tech- 
nique well suited for the semiclassical domain we are interested in, to be rather useful in this respect. In fact, 
many years ago a wave-mechanical description of the collisionless ion dynamics around electric probes has been 
pursued by Liu ll35l but he found no followers. 



4 Charging model 



Using the results of the previous section, a model for the grain charge taking surface states into account can be 
constructed as follows. Within the disturbed region of the particle, the density of free electrons (ions) is much 
smaller than the density of bound electrons (ions). In that region, the quasi-stationary charge (again in units of 
— e) is thus approximately given by 



Z{x) = AtiR^ I dx'{l+x') 



nl{x') - n\{x') 



(11) 



with X < \f = y^kTi/Aire'^ni, the ion Debye length, which we take as an upper cut-off, and ^ the density 
of bound electrons and ions. For the plasma parameters used in Fig. 121 Af ~ 60fim. The results for the surface 
states presented above suggest to express the density of bound electrons by an electron surface density: 



n^{x) ~ aeS{x — Xe)/R 



(12) 



with Xe — Xb — and the quasi-stationary solution of Eq. ^ without the recombination term. Equation ^ 
is thus still a rate equation on the grain surface. We will argue below that once the grain has collected some 
negative charge, not necessarily the quasi-stationary one, there is a critical ion orbit at a;, 1 — 10 ^ Xe 
which prevents ions from hitting the particle surface. Thus, the particle charge obtained from Eq. STT[ is simply 
Zp = Z{xe < X < Xi). Inserting Eq. (ITSl i into Eq. (fTTl i and integrating up to x with Xe < x < Xi leads to 
Eq. (|5]l, the expression for the particle charge deduced from the rate equations (|2]i and (jS} under the assumption 
that ions do not reach the grain surface on the microscopic scale. 

For an electron to get stuck at (to desorb from) a surface it has to loose (gain) energy at (from) the surface ll28l . 
This can only occur through inelastic scattering with the elementary excitations of the grain schematically shown 
in the left panel of Fig.|3l To calculate the product (sr)p requires therefore a microscopic model for the electron- 
grain interaction. First steps in this direction were taken in ll30l . 

In Ref. Il25l we invoked however the phenomenology of reaction rate theory to approximate (sr)e. Specif- 
ically, we assumed that electrons with rather low and rather high energies are, respectively, reflected by the 
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grain boundaiy grain boundary 

Fig. 3 Left panel: Schematic illustration of electron energy relaxation at the grain boundary. The electron looses (gains) 
energy due to creation (annihilation) of elementary excitations of the grain. Due to these processes it may get trapped in 
(escape from) the discrete bound states of the short-range polarization potential. Right panel: Schematic illustration of ion 
energy relaxation in the vicinity of the grain surface. The ion looses (gains) energy due to collisions with other plasma 
particles (electrons, ions, neutrals). It is due to these processes that an ion may get trapped in (escape from) the continuum of 
bound states of the long-range Coulomb potential. 



Coulomb and surface barrier of the particle and that sticking (desorption) primarily affects electrons with ener- 
gies slightly above U. After overcoming the Coulomb barrier this group of electrons is almost in equilibrium 
with the surface electrons. We can thus apply the Lennard-Jones-Devonshire formula 12811361 to obtain, 



h 



knT, 



exp 



E: 



•d -\ 



(13) 



where h is Planck's constant, Tp is the grain temperature, and is the electron desorption energy, that is, the 
binding energy of the surface state from which desorption most likely occurs. The great virtue of this equation 
is that it relates a combination of kinetic coefficients, which depend on the details of the inelastic (dynamic) 
interaction, to an energy, which can be deduced from the static interaction alone. Kinetic considerations are thus 
reduced to a minimum. They are only required to identify the relevant temperature and the state from which 
desorption most probably occurs. 

Equation (|5]l is a self-consistency equation for Zp. Combined with Eq. ( fTsl l. and approximating the electron 
plasma flux jpiasma orbital motion limited flux. 



27rm, 



■ exp 



~Zpe^ 



RksTe 



(14) 



which is reasonable, because, on the plasma scale, electrons are repelled from the grain surface, the grain charge 
is given by 



ksTp 



exp 



E^ 



ksTp 



J^^Zp) . 



(15) 



Thus, in addition to the plasma parameters rie and Te, the charge depends on the surface parameters Tp and E'^. 

Without a microscopic theory for the inelastic electron-grain interaction, a plausible estimate for Ef has to be 
found from physical considerations alone. Since by necessity the electron comes very close to the grain surface it 
will strongly couple to elementary excitations of the grain (see Fig|2]and the left panel of Fig. |3]l. Depending on 
the material these may be bulk or surface phonons, bulk or surface plasmons, or internal electron-hole pairs. For 
any realistic surface barrier, where the electron wavefunction leaks into the solid, the electron will thus quickly 
relax to the lowest surface bound state B The n — 1 state for the infinitely high barrier is an approximation to 



The microscopic model for electron energy relaxation at metallic boundaries employed in (30| works also for an infinitely high barrier. 
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Table 1 Plasma parameters 12411331 used to obtain the results plotted in Fig. [4] 



P[W] p[Pa] ne,n,[109cm-3] kTejeV] kTi[eV] 



5 


102 


0.33 


1.5 


0.030 


12 


102 


0.62 


2.2 


0.040 


20 


102 


0.9 


1.7 


0.035 


30 


102 


1.2 


1.28 


0.036 


40 


102 


1.5 


1.2 


0.038 


50 


102 


1.4 


1.0 


0.039 


60 


102 


1.6 


0.8 


0.040 


12 


30 


0.26 


2.2 


0.030 


12 


40 


0.4 


2.2 


0.030 


12 


60 


0.48 


2.2 


0.030 


12 


80 


0.52 


2.2 


0.036 


12 


102 


0.62 


2.2 


0.040 


12 


120 


0.8 


2.2 


0.040 



that state. Thus, it is reasonable to expect 

where is the lowest eigenvalue of the electronic Schrodinger equation. For an MF particle with e = 8 this 
leads to ~ O.beV. The particle temperature cannot be determined in a simple way. It depends on the balance 
of heating and cooling fluxes to-and-fro the particle and thus on additional surface parameters ll37l . We use Tp 
therefore as an adjustable parameter. To reproduce, for instance, with Eq. (fT5] l the charge of the particle in Fig.|2] 
Tp = 395 K implying (sr)^ ~ IQ-^ s. 

Equation ( fTSl l depends on the assumption that once the particle is negatively charged ions are trapped far away 
from the grain surface. Indeed, treating trapping of ions in the field of the grain as a physisorption process suggests 
this assumption, which is perhaps counter-intuitive. Provided the ion remains intact on its way towards the grain 
surface (no ion neutralization in the disturbed zone of the grain due to electron capture) it gets bound to the grain 
only when it looses energy. Because of its low energy and the long-range attractive ion-grain interaction, the ion 
will be initially bound very close to the ion ionization threshold (see Fig. |2]i. The coupling to the elementary 
excitations of the grain is thus negligible and only collisions with other plasma particles are able to push the 
ion from an extended state with positive energy to a bound state with negative energy and then from a given 
bound state to a lower one (see the illustration in the right panel of Fig. [3]i. Since the interaction is classical, 
collisions, for instance, charge-exchange scattering between ions and atoms, act like a random force. Ion energy 
relaxation can thus be envisaged as a destabilization of orbits. This is in accordance to what Lampe and coworkers 
assume [20 2A^22i- In contrast to them, however, we ll25l expect orbits whose spatial extension is smaller than 
the scattering length to be stable because the collision probability during one revolution becomes vanishingly 
small. For a circular orbit, a rough estimate for the critical radius is ri~R{\ + Xi)~ (27r(Tcxfig)~^ which leads 
to a^i ~ 5 ^ Xe ~ when we use the parameters of the helium discharge of Fig.|2]and CTcx = 0.32 x 10~^^ cm? 
which is the measured cross section at 0.3 eV ll38l[39l . 

Indeed, Lampe and coworkers approach 1201 1211 l22l shows a pile-up of trapped ions in a shell of a few 
radius enclosing the grain. They would however not expect a relaxation bottleneck. This point can be only 
clarified with a detailed investigation of the ion dynamics and kinetics in the disturbed region of the grain tak- 
ing the complete kinematics of charge-exchange collisions encoded in the differential collision cross section and 
the centrifugal barriers separating bound from unbound ion motion into account. In fact, Lampe and coworkers 
neglect the momentum transfer during a charge-exchange coIUsion as well as the barriers. Tskhakaya and cowork- 
ers 1261 l27l on the other hand pointed out that the latter could severely overestimate the collision-enhanced ion 
flux. In reality, this flux, they claim, is much smaller than the one obtained by Lampe and coworkers. If this is 
indeed the case, the charges obtained from the collision-enhanced ion flux model would be much closer to the 
orbital-motion limited ones and thus far away from the experimentally measured charges (see next section). 
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Fig. 4 Left and right panel show, respectively, the power and pressure dependence of the charge of a MF particle {R = 
4.7 nm) in the helium discharge of Ref. 1331 . The particle temperatures reproducing the experimental data (filled squares) 
are indicated and the charges obtained from j^^^ = jf (OML) and = + jp^ (OML & CX) with tJcx = 

0.3 X lO"^** cm 1381 1391 . which is the measured cross section at 0.3 eV, are also shown. The predicted increase of the 
particle temperature with pressure is plotted in the inset of the right panel. 



In ll25l we pushed the assumption of a critical ion orbit to its Umit and approximated the density of ions 
accumulating in the disturbed region of the grain, and being responsible for the partial screening of the grain 
charge, by a surface density cr; which balances at Xi the ion charging flux with the ion desorption flux (see right 
panel of Fig.[T]i. Mathematically, this gives rise to a rate equation similar to (O, but without the recombination 
term and interpreted as a rate equation at r = r^. At quasi-stationarity, the ion surface density is thus (T; = {sT)iji. 
Although Eq. ( fT3b assumes elementary excitations of the grain to be responsible for sticking and desorption we 
expect a similar expression (with E^, Tp replaced by Ef, Tg) to control the density of trapped ions. Equation (fTTI) 
leads then to Z{xi < x < \f) ^ Zp - Zi with 



=4^i?2(l+X,)2— -CXp 

ksTg 



EfiZp) 



jf (17) 



ksTg 

the number of trapped ions where we assumed that the critical ion orbit, which is near the sheath-plasma boundary, 
is fed by the Bohm ion flux jf^ = .Qni^J ksTe / nii. 

The ion desorption energy is the negative of the binding energy of the critical orbit, 

Ef(Zp) = -V,{x,)U = ATTdc^aBngZpRo , (18) 

and depends strongly on Zp and Xi. For the situation shown in Fig.|2]we obtained, for instance, using Tg = Tp = 
395 K, the particle temperature which reproduces Zp = 6800, Ef ~ 0.37 eV and (sr)j ~ 0.6 x 10"** s. The 
ion screening charge is then ~ 148 <C Zp which is the order of magnitude expected from molecular dynamics 
simulations BOl . Thus, even when the particle charge is defined by Z{xi < x < Xf) it is basically given by Zp. 

5 Results 

We now use Eq. ( fTSl ) to calculate for various experimental set-ups the charge of a dust particle. In all cases the 
plasma parameters are known. The only free parameter is thus the particle temperature. Although there exist 
optical methods to measure the temperature of dust particles BTl they have not yet been used in conjunction with 
charge measurements. In view of the importance of the particle temperature for the surface model, we hope that 
in the near future an experimental group can be found to perform such measurements. 

In addition to the results obtained from the surface model we also plot charges deduced from balancing the 
orbital motion limited electron flux j^^^ given in Eq. (fT4l) with, respectively, the orbital motion limited ion flux. 



J 



OML ' knT 



27rmj 



Zpc^ 



RksT, 



(19) 
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Fig. 5 Left panel: Pressure dependence of the charge of a MF particle with R = 1 in the neon discharge of Ref. II12I 
(squares) 1251 . Middle and right panel: Radius dependence of the charge of the MF particle for p = 50 Pa and p — 100 Pa, 
respectively. In all panels solid lines denote the (isothermal) charges deduced from the surface model whereas dotted, dashed, 
and long-dashed lines are the charges obtained from balancing on the grain surface j^^^ with ji'^^+ji''^ using, respectively, 
o"cx = 0.41 X lO^^^cm^ 1421 . which is the experimentally measured cross section at 0.12 eV, (Tcx ~ 1-0 x 10" 
CTcx = 2.0 X 10""cm^. 



and the collision-enhanced ion flux (see ||201|2T1 and also lfT2l ). 

o2 \D 



•OML , -ex ^ I ksTi 



RkBTi ' lex \RkBTi I 



(20) 



where /cx = (ccx'^g)^^ is the scattering length, ctcx is the charge-exchange cross section, and rig = p/ksTg is 
the gas density. 

We start with MF particles confined in a helium discharge. In Fig. |4]we show the power and pressure de- 
pendence of the charge of a particle with i? = 4.7 fim at rest in the helium discharge of Ref. [|33l, the plasma 
parameters of which are given in Table [T] The p = 102 Pa data point of the P = 12 W run served as an illus- 
tration in Fig. ID Using the parameters of Table [T] we calculated from Eq. ( fTSI ) for each data point Tp such that 
Zp{Tp) = ^oxp- The power dependence of the charges, shown in the left panel of Fig. H] could be reproduced 
by a single particle temperature Tp = 405 K while the pressure dependence shown in the right panel required 
to adjust for each pressure the particle temperature. Assuming Tp to scale with the gas temperature Tg, the pre- 
dicted particle temperature plotted in the inset of the right panel is in accordance with what we would expect 
from the pressure dependence of Tg in noble gas discharge ifjTl indicating that our approach gives physically 
consistent results. For comparison we also plot the particle charges deduced, respectively, from j^^^ ~ Ji'^^ 
and j^^^ = jP^^^ + iF"^. Obviously, the agreement with the data is not very good. 

That the charges obtained from the orbital motion limited flux balance are not too close to the experimental 
data is expected. But even the charge-exchange enhanced flux balance gives not particularly good results. This 
can be also seen in Fig.|5]where we analyze the charges of MF particles confined in the bulk of the neon discharge 

of Ref. na. 

The pressure dependence of the charge of a MF particle with R = 1 /.im is shown in the left panel of Fig.|5] 
Since the plasma parameters entering Eq. ( fTSl ) are known il21 . Tp is again the only free parameter Fixing Tp at 
a particular value gives the isothermal particle charges Zp{Tp) shown by the solid lines. From Zp{Tp) = Z^xp 
follows then the Tp required to reproduce the experimental charge. The predicted increase of Tp with pressure is 
again in accordance with the results of the calorimetric study of noble gas discharges presented in ll37l . 

In Fig. 12] we also plot the charges obtained from the charge-exchange enhanced flux balance condition. For 
CTcx = 2 X 10"^'' cm^ the agreement with the data is in fact quite good but this value of the cross section is 
almost five times larger than the experimentally measured value (ctcx = 0.41 x 10^^^ cm^ at 0.12 eV 1421 ). 
That the agreement with the experimental data is not too good can be also inferred from the radius dependence 
of Zp for p = 50 Pa and p = 100 Pa which we show, respectively, in the middle and the right panel of Fig. |5] 
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Fig. 6 Radius dependence of the charge of a MF particle in the sheath of an argon discharge at p = 6.67 Pa (left panel) 1251 
and p — 13.34 Pa (right panel). Squares are experimental data from 1101 and solid and dashed lines give, respectively, the 
charges obtained when the depletion of Ue in the confining sheath is included or not. For comparison we also show the charges 
deduced from = jf^^ (OML) and from = + jf^ (OML & CX) with cTcx = 0.72 x IQ-^^cm? which 

is the experimentally measured cross section at 0.1 eV 1421 . 



Clearly, the radius dependence of the grain charge seems to be closer to the nonlinear dependence obtained from 
Eq. ( fTSl l than to the linear dependence resulting from the charge-exchange enhanced flux balance, irrespective of 
the chosen value of the charge-exchange cross section. 

Finally, Fig. |6] showing for two different pressures the radius dependence of Zp for MF particles confined 
in the sheath of an argon discharge fTOl, the plasma parameters of which are given in table |2] provides addi- 
tional support for our model. To approximately account for the fact that particles with different radius expe- 
rience different plasma environments, we included the depletion of Ue in the sheath by replacing rie in j'^^^ 
by rie exp[e$(zcq(-R))/A:BJe] with $(2) the sheath potential and Zcq(^) the measured equilibrium position of 
the particle with radius R lITOl . When the grains are not too deep in the sheath (R < 5 fim), we find for 
p = 6.67 Pa and p = 13.34 Pa excellent agreement with the data for, respectively, Tp ~ 410 K — 430 K 
and Tp = 400 K — 420 K, although the particle temperatures are perhaps somewhat too high in view of the 
small amplitude of the rf voltage. Our approach fails, however, completely for R > 5 /im. We attribute this to 
the ad-hoc modification of the electron flux which may not capture the total electron flux close to the electrode. 
An improved treatment would calculate the electron flux self-consistently taking not only the depletion of the 
electron density into account but also the flux due to sub-thermal secondary electrons from the electrode 1111 . 



6 Conclusions 

In this paper, we discussed the main assumptions underlying the surface model proposed in ||25l for the calcu- 
lation of the charge of a dust particle immersed in a quiescent plasma and confronted the model with additional 
experimental data. 

The main hypothesis of the model, suggested by the analogy of charging of dust particles with physisorption of 
charged particles to the particles' surface, is that electrons and ions are, on a microscopic scale, spatially separated 
because the potential in which physisorption of electrons takes place is the attractive short-range polarization 
potential whereas for ions it is the attractive long-range Coulomb potential. For electrons the analogy may be 

Table 2 Plasma parameters ||9) used to obtain the results presented in Fig.|6] 

p[Pa] nejlO^cm-^] kTe[eV] m[eV] Vri\V] 
6.6 1.7 3.7 0.026 96.4 

13.6 2.4 3.9 0.026 96.4 
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obvious. But for ions it may be not. However, without the grain no attractive long-range Coulomb potential would 
exist. Thus, the ion dynamics and kinetics in the vicinity of the grain is also a kind of surface physics although it 
occurs deep in the disturbed region of the grain and may be strongly affected by the plasma environment. 

Within the surface model, the grain charge and its partial screening can be calculated by balancing, on two 
different surfaces, the electron collection flux with the electron desorption flux and the ion collection flux with 
the ion desorption flux. The charge of the grain is then given by the number of electrons quasi-bound in the short- 
range polarization potential of the grain whereas its partial screening is given by the number of ions quasi-bound 
in the long-range Coulomb potential. 

The grain temperature turns out to be an important parameter Since, for the experiments we analyzed, the 
grain temperature has not been measured we used it as an adjustable parameter and obtained for physically mean- 
ingful particle temperatures excellent agreement with the experimentally measured grain charges. We challenge 
therefore experimentalists to simultaneously measure the grain charge and the grain temperature. Our model 
could then be easily tested. 

The charges obtained from the surface model depend on the surface parameter (sT)e which is the product of 
the electron sticking coefficient with the electron desorption time. Both parameters depend on the inelastic, that 
is, the dynamic interaction between the electron and the grain. As discussed in 1301 . a rigorous calculation of 
(sT)e has to be based on a microscopic model for the electron-grain interaction taking elementary excitations 
of the grain into account. Quantum-kinetic equations have then to be solved to obtain Sg and Te separately. 
For the product, however, a rough estimate, which turns out to be surprisingly good, can be obtained from the 
Lennard-Jones-Devonshire formula (fT3T l relating (sr)e to the particle temperature and the electron desorption 
energy. Measuring these two quantities directly would eliminate any free parameter from our model. 

The surface model is a first attempt to treat plasma-controlled electron and ion fluxes and material-controlled 
plasma-wall interactions at the grain surface on an equal footing. Even in the present rudimentary form the 
model performs better than approaches relying exclusively on an improvement of the plasma fluxes. Hence, if 
nothing else, it indicates that the charge a dust particle acquires in a plasma depends not only on the macroscopic 
plasma environment but also on microscopic processes on the surface occurring on a scale which is beyond the 
Boltzmann-Poisson description of the plasma-grain interaction. A quantitative theory of grain charging has to be 
therefore based on a systematic exploration of this ultimate boundary layer. 

Acknowledgements Support from the SFB-TR 24 "Fundamentals of Complex Plasmas" is greatly acknowledged. F. X. B. 
is also grateful to M. Lampe for a particularly illuminating discussion. 
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